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Abstract 


In  this 

paper 

we . 

discuss  the 

quantum  coherence 

properties  of 

the 

FEL 

within  the 

f  r  a  m  ework  of  the 

single  electron  n o n- r e 1  a t i v i s t  i  c  Hamiltonian  picture. 
We  analyse  the  problem  both  in  the  single  and  multi¬ 
mode  hypotheses  . 


( + )  Quantum  Institute,  University  of  California  at  Santa 

Barbara 

(++)  ENEA  Student 

■'-M-e)  T.N.F.N.  Sez.  Napoli 
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The  quantum  non-relativistic  Hamiltonian  picture 

[1]  of  the  Free  Electron  Laser  (FEL)  has  provided 

(*) 

a  useful  tool  to  understand  the  coherence  proper¬ 

ties  of  this  new  kind  of  laser  device  [2-5]. 

In  this  paper  we  will  review  the  main  features 
of  that  analysis,  relevant  to  a  single  mode  -  single 
electron  picture.  Furthermore  we  show  how  the  longi¬ 
tudinal  multimode  effects  can  be  properly  accounted 
for,  in  a  quantum  mechanical  way. 

The  basic  elements  of  the  non-relati vistic-quan- 
tura  Hamiltonian  Theory  of  the  FEL  will  not  be  reported 
here,  but  they  can  be  found  in  many  review  articles 
(as  e.g.  Ref.  [6]  where  a  detailed  account  of  both 
classical  and  quantum  FEL  treatments  has  been  given). 

The  non-relativistic  single  mode  -  single  elec¬ 
tron  FEL  Hamiltonian  reads 


+  ilu)(aLaL+l/2) 


+  Tiw(ayau+l/2) 


_or  +  -2ikz  +  2ikz 
+  T\«[aLaue  +  auaLe 


(1) 


(For  the  symbols  see  Table  I) 


The  physical  meaning 

of 

( 1 )  is 

transparent.  The 

first  three  terms  in 

(1) 

account 

for 

the  electron 

energy  and  free  field 

ones ,  while 

the 

FEL  mechanism 

arises  as  due  to  the  fourth  term  in  the  Hamiltonian, 
which  causes: 


(*)  From  now  on  when  we  say  coherence  we  refer  to  quantum  co¬ 
herence  and  not  to  the  spatial  and  temporal  ones. 


a)  The  creation  of  a  laser  photon,  the  destruction 
of  an  undulator  one  and  the  loss  of  2TiK  of  electron 


momentum . 

b)  The  conjugate  process,  namely  the  destruction 
of  a  laser  photon,  the  creation  of  an  undulator 
one  and  the  gain  of  2HK  of  electron  momentum. 

It  is  well-known  that  [6]  the  Hamiltonian  (1) 
allows  two  laws  of  conservation  relevant  to  the 
total  number  of  photons  (laser  +  undulator)  and 
to  the  total  momentum  (electron  +  fields). 

The  above  denoted  conservation  laws  suggest 
that  the  evolution  of  the  FEL  operation  states  can 
be  characterized  by  a  single  integer  only,  namely 
the  number  of  exchanged  photons  during  the  inter¬ 
action  . 

In  other  words,  if  we  start  from  a  single  state 
with  electron  energy  P2  /2m  and  n°  ^  initial  photon 
number,  the  following  states  only  are  coupled 


|4>>  = 

exp 

p2 

“K  +  nu  * 

i  >  1 1  j  • 

• 

ao 

E 

O 

rt 

Q 

-  22.Hk ,  n°L 

+  nj- 

£  >  (2) 

*<=> 

ii 

o 

where 

l 

is  the 

number  of 

exchanged 

photons  and 

Cz(t)  are 

time  dependent  co e f f ici en t s . 

If 

we  assume 

that  the 

initial  ’’ 

laser"  field 

is  the 

vacuum  (n° 

L 

=  0)  ,  we 

find  that 

the  coef  f i- 

cients 

obey  the 

f ol lowi ng 

differential-difference 

equation 


4 

i  dC^  =  (-wQ  +  e£  H'C  +  TT[  /(5.  +  1  )  ( nj’j-S,  )  + 

d  T 

♦  /l(nS  -l+1)  VJ  •  V0)  ■«!,.  (3) 

The  above  expression  is  known  as  the  Spherical  or 
SU2  -  Raman-Nath  (R.N.)  equation  and  has  been  recently- 
discussed  in  Refs  [7].  It  cannot  be  solved  exactly 
by  means  of  known  functions,  but  a  perturbed  analysis 
in  terms  of  the  parameter  e ,  linked  to  the  quantum 
electron  recoil  (see  Table  I),  can  be  accomplished. 
Let  us  notice  that  at  the  zero*"^  order  in  e  ,  and  thus 
in  the  quantum  corrections,  the  solution  of  (1)  can 
be  written  quite  straightforwardly  [7],  so  that  the 
probability  of  emitting  l  photons  in  the  recoilless 
approximation  is  given  by 

|C° | 2  =  (  U)  P*(l  -  P)(nU~Z)  (4) 

where 


5  =  v^w 

l  +  * 

(5) 

It  is 

easy  to  recognize  (4)  as  a 

binomial  distribu- 

t  ion . 

However,  since  the  number 

of  undulator 

"pho- 

tons" 

is  always  a  large  quantity, 

taking  the 

limit 

n°>H 

(4)  gives 

2 
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i  £  =0,  n  °>>Z 


--  £  !3(t)i^e'|a(T)l 


flt(t)  =  (-i)expj  iw0x/2 


a 


sin( w0 t/2) 
wn/2 


(6) 


Needless  to  say,  the  above  expression  (6)  is  a  Poisson 


distribute  on  , 

therefore 

up 

,  th 

to  the  zero 

order  in 

the  quantum 

corrections 

we 

have  recovered 

a 

result 

strongly  reminiscent  of 

the 

conventional  laser 

stati- 

sties. 

We  could 

also  show, 

in 

a  more  rigorous 

way , 

wi  th- 

in  the  framework  of  the  same  approximation,  that  the 

states  (2)  evolve  from  the  vacuum  into  coherent  G 1  a  u  - 
(*) 

ber  states 

It  is  easy  to  realize  that  the  recoilless  approxi¬ 
mation  amounts  to  a  gainless  process  [4,8].  The 
inclusion  of  e  is,  therefore,  necessary  if  we  want 
to  find  meaningful  results.  If  we  look  only  for 
first  order  quantum  corrections  we  derive,  however, 
significant  deviation  from  the  previously  recovered 
results  [2,4,8].  Namely 


a)  We  cannot  prove  that  the  states  (2)  evolve  from 
the  vacuum  into  Glauber  ones. 

b)  The  probability  of  emitting  l  photons  is  no  longer 
a  Poisson  situation  but  is  given  by 


I2 

£  *  0 


'ir 


e-k(T)|2|a(T)|U 
Z  ! 


1  +  |  ot  (t  )  1 


|  a  (x  )  |  •  [(2£  +  l )  |  a  (t  )[2-  a  (t  )  |  4  -£2  ]  +  .  .  . 


(*) 


For  a  more  detailed  account  see  Refs  [4], 
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c)  Furthermore  we  can  expect  genuine  quantum  effects 
such  as  antibunching  and  squeezing. 

Evaluating,  indeed,  the  quantity 

<A£2>-  <l>  =  -  e  ~~  |  Qt  ( T  )  ]  **  (8) 

3  w0 

We  find  that,  in  the  FEL  gainless  process,  the  photon 
number  fluctuations  are  just  the  square  root  of  the 
number  of  the  emitted  photons,  as  in  conventional 
,  while,  if  there  is  gain  or  loss,  bunching 
or  antibunching  arise.  (See  also  Refs  [2,4,8]  for 
further  .comments.) 

As  to  the  squeezing,  if  we  introduce  the  operators 
A,  »  ~(A  +  +  A)  (A+  =  a  +  e"2ikZ) 

A2  -  7~(A+  '  A)  (9) 

which  realize  a  kind  of  electron-field  p,q  variables 
we  can  evaluate 

<AA  l>  =  j  +  ~  z  ( ~  |  ot(  T )  |  2  )cosw0T  +  ... 

4  z  3  w0 

11  9 

<AAj>  --  -  ~  e (- - |a(x) | 2 )cosw0x  +  ...  (10) 

4  Z  Cf  W  q 

We  must,  therefore,  underline  that  both  antibunching 
and  squeezing  are  related  to  the  gain  process  and 
are  the  only  genuine  quantum  effects  which  can  be 
derived  within  this  analysis. 
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The  above  results  hold  in  the  hypothesis  that 
the  initial  laser  field  is  the  vacuum. 

Let  us  discuss  the  case  of  an  arbitrary  coherent 
Glauber  state  as  initial  "laser"  field.  This  is 
the  physical  situation  of  an  amplification  FEL  experi¬ 
ment  and  we  are  going  to  discuss  in  which  way  the 
coherence  properties  of  a  conventional  laser,  under¬ 
going  amplification,  may  be  affected  by  the  interac¬ 
tion. 

To  carry  out  this  analysis  we  need  only  a  few 
modifications  to  the  above  procedure.  Since  we  do 

not  start  from  a  state  with  a  fixed  number  of  photons, 

( *) 

but  from  a  coherent  one  we  must  rewrite  (2)  as 
l  P 2 

j  =  expj-i[— 2— +  u>(n°+l)]t 

l  2mtl  u 


00  CO 

•  L  L  er'  — |2/2<*o  cj(t)|n;  P0-  l> 

n=o  £  =-  a  ^n  ! 

In  this  case  the  differential  equation  for 
coefficients  writes 


(ID 

the  C. 
£ 


ip11  __ 

i  —l  =  (-w0+  e£  )£c"  +  8  [/n  +  £  +  1  c"  + 

dx  £  R  £  +  i 

w77T0'  c"<0)  ■  V.  (12) 

The  above  expression  is  known  as  the  Harmonic  R.N. 
equation  recently  analysed  in  Refs  [9].  Even  in  this 


(*) 


From  now  on  we  assume  that  the  number  of  undulator  "photons" 
be  very  large  and  that  its  variation  can  be  neglected,  fur¬ 
ther  comment  on  this  assumption  will  be  given  in  the  fol¬ 
lowing. 
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case,  an  exact  solution  cannot  be  found  in  terms  of 
known  functions;  but,  as  before,  we  will  carry  out 
only  a  first  order  analysis  in  terms  of  the  electron 
recoil . 

The  zero*"^  order  solution  can  be  easily  found 
(9)  and  writes 

fT 

C?  =  exp  i— !  ot  (  T  ’  )  |  2  dr'  <J>  [a(x)]  (13) 

x#  o  2  n 

o 

where 


‘ /(7+ijT  ='|a<T)|!/2(“(^)i^[|a(t)| 


with  L^(*)  generalized  Laguerre  polynomials. 

We  must  now  stress  a  significant  difference  with 
respect  to  the  previous  case.  Even  in  the  recoilless 
approximation  we  cannot  prove  that  (11)  is  a  Glauber 
state  at  any  time. 

This  is  due  to  the  fact  that  the  creation  annihi¬ 
lation  operators,  acting  on  our  Hilbert  space,  charac¬ 
terized  by  the  integer  2,  only,  affect  both  the  fields 

(*) 

and  electron  variables  see  (9), (-8),  so  that 


(a(r) 


+  ane 


2ikz 


(14) 


The  relationship  (14)  does  not  define  a  coherent 
state,  because  is  not  a  c-function  but  an  opera- 


(*) 


If  A  were  a  pure  field  operator  we  could  easily  prove 


A|^>  =  (a(t)  +  a0)|i|» 

£=0  “  £=  0 

and  thus  Glauber  coherence 
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tor.  Notwithstanding,  we  can  still  prove  that,  at 
the  zero'"'1  order  in  £  ,  we  can  recover  the  well-known 
results  of  the  Poisson  statistics.  However,  if  we 
want  to  recover  the  gain  we  must  look  for  a  first 
order  perturbed  expansion  in  the  recoil.  Therefore, 
including  the  lowest  order  quantum  corrections,  the 
evolution  of  the  states  (11)  reads 


exp 


|  [  a(t  '  )|  2  dt  ’ 


+  icuAt(n°  +  l)T 


+ 


.  Po 

+  i — rAtt 
2mti 


E  •' 


2  /  2  n 
<*o 


1 

(  n  +  £)  ! 


[  a(  t )  1  ^  • 


-  I  a(  x)  I  2 / 2 


e 

where  A .  and  D  . 

2 ,  n  x.,  n 

sions  given  by 


[A,  +  iD.  ]  I  n  .  P  -22tik  ;  2>  (15) 

are  somewhat  complicated  expres- 


.  x  e  i sXil 

2 ,  n  n  j  ct  (  T  )  1  3w0 


(22  +  l)L^  +  1 ( • ) ! a ( T  )  |  2 

n 


-  (22-l)(n+2)L^_1  (  *  ) 
n 


e  3  1  a(r  )  1 
a  (  j )  |  3  w  0 


a  ,  2  —  2 

|  at  (  T  )  |  4  L  (•)  -  (n+2)(n+2-l)L  (•) 

n  n 


D 


2 ,  n 


£  T 
2 


£  +  1  £  —  1 

( 22+1 ) | a(x ) I  2  L  (*)  +  (  22  - 1 ) ( n+2  )  L  (.) 

n  n 


-  |  a(  T)  |  2(  • )  +  (  n  +  2)  (  n+2 -1  )  [}  *(•) 

2  n  n 


(16) 


The  above  quantities  can  be  exploited  to  evaluate 
a  number  of  interesting  quantities  such  as  the  average 
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photon  number,  the  second  moment  of  the  distribution 
and  so  on.  The  explicit  calculation  of  these  quanti¬ 
ties  is  tedious  but  straightforward  it  requires  a 
repeated  application  of  the  well-known  properties 
of  the  Laguerre  polynomials  [10];  so  that  after  some 
algebra  we  find 

<n+1>£?<o  =  ^a°'2  +  ~  £  t 1  +  2  1  a  0  ]  2  ] 

(17) 

It  is  easy  to  recognize  the  physical  meaning  of  the 
various  terms  appearing  in  (17).  The  first  terra 
is  the  average  number  of  photons  of  the  input  field, 
the  second  one  accounts  for  the  spontaneously  emitted 
photons  and  finally  the  third  term  is  the  stimulated 
part  (the  extra  term  1  comes  from  the  vacuum  field 
fluctuations)  . 

We  can  also  evaluate  the  quantity 


<  A  (  n  +£  ) 2  >  ,  -  <n  +  £>  =  -e  (2|a0|2+l) 

E  *  o  s*  o  101  3  w  . 


-  2  e | a0 


2 


3  |  a  1  2 
3  w0 


(18) 


Thus  finding  again  that  we  can  have  bunching  or  anti¬ 
bunching  according  to  the  case  w0  <  0  (see  Fig.  1). 
Similar  conclusions  can  be  drawn  for  the  squeezing. 

Before  discussing  the  problem  of  the  multimode 
analysis  we  need  some  comments  on  the  assumption  of 
constant  undulator  photon  numbers  during  the  interac¬ 
tion.  It  is  well-known  that  the  concept  of  undulator 
,  ho  ton  is  grounded  on  the  Weizsi'cker -Williams  approxi- 
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mation  [6].  Furthermore  it  is  also  well-known  that 

an  undulator  with  a  static  field  of  1  kG  has  an  equi¬ 
valent  photon  number  density  of  about  102I*/cm3,  which 

is  so  large  that  it  is  practically  unaffected  by  the 
interaction.  Our  assumption  is  therefore  well  justi¬ 
fied.  Notwithstanding  it  is  not  correct  for  a  proper 
and  rigorous  analysis  of  the  FEL  coherence.  We 

should  indeed  treat  the  undulator  field  as  a  coherent 

Glauber  one,  in  this  way  we  reduce  the  problem  of 
FEL  coherence  to  that  of  the  evolution  of  two  coupled 
coherent  fields.  This  problem  may  be  understood 
in  terras  of  the  evolution  of  the  so-called  coherent 
Bloch  states  [12].  In  this  case  even  when  we  start 

from  a  zero  laser  field  we  meet  with  the  problem  of 

the  electron-field  operators,  discussed  before.  There¬ 
fore,  strictly  speaking,  the  problem  of  the  FEL  Glau¬ 
ber  coherence  even  at  the  zero^  order  in  e  is  still 

open. 

The  main  implicit  assumption  of  the  up  to  now 

developed  analysis  is  that  the  electron  has  a  fixed 

initial  momentum  and  that  it  is  completely  delocalized. 

We  will  now  assume  that  the  initial  electron  wave- 

-function  is  an  arbitrary  packet  and  that  the  laser 

field  may  be  expanded  into  many  longitudinal  modes 

(*) 

[12].  Therefore  we  have  that  the  Hamiltonian  now 

writes 

p  2  n  ^ 

H  =  —  +  hu)..(a  +  a  +1/2)  +  ft  V"  u>  (a  +  a  +1/2)  + 

2m  UUU  * —  sss 

s=  1 

(*) 

A  proper  analysis  of  the  FEL  signal  starting  from  the  vacuum 
requires  a  multimode  expansion  even  in  the  hypothesis  of 
continuous  electron  beam. 
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»  E 


a 


S,  3 
j<s 


a+a  exp [ i ( k  -k.)z] 
s  J  s  j 


+  h  .  c  . 


(19) 


The  physical  meaning  of  (19)  is  self-explanatory. 

We  have  a  further  terra  namely  the  fourth  one 
which  accounts  for  the  laser-laser  interaction  (for 
further  comments  see  Ref.  [12]). 

According  to  the  Hamiltonian  (19),  to  the  assump¬ 
tion  of  arbitrary  wave  packet  for  the  electron  and 

(*) 

to  the  laws  of  conservation  we  can  write  the 

analogous  of  the  states  (2)  as 


u(o>  =  Ej  «*<V<Wvt)  • 

1  &  J 


•  exp 


j-  xfi  At  +  n°u)  At  +  r  .At]t  I  • 
1  2m  U  U  J  J  1 


n  u 

|  k  -  Y'  (k  +k  )l  ,  [n°+2,  }  ,  n°  -  i  > 

e  s  u  s  s  s  U  t—  s 


s  =  i 


s=  1 


(20) 


It  is  also  easy  to  prove  that  the  coefficient 
C  (k  ,  t)  may  be  determined  from  the  following 
equations 


dCfk^r) 


dt 


=  Y  [-  q  +  l  e  l  12.  .C,  .  (k  ,  x )  + 

J  s  s  ,  j  s  j  lx.)  e 


+  Y  (n .  ,. ) 

j-  j'u 


(n°+i  .  +  1  )  (  n  0  -  LI  )  C  ,+ 

J  J  U  q  q  U  .  ) 

J+i 


(*) 


The  laws  of  conservation  are  again  relevant  to  the  total 
momentum  and  total  number  of  photons.  See  Ref.  [12]  for 
the  specific  details. 
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The  above  equation  is  an  example  of  SU  -R.N.  [14] 

n 

equation  we  do  not  dwell  (20)  on  its  general  form 
but  we  make  some  simplif icative  assumptions. 

(1)  We  neglect  the  undulator  photon  number  variation 
and  start  from  the  vacuum. 

(2)  We  neglect  the  laser-laser  interaction. 

So  that  (20)  reduces  to 


dC 


{£} 


d  x 


V 

r  [  -  n .  +  e  .111.  c 

r~,  J  s  s,j  s  j  {£} 


i=  i 


E  nRj  \-/lj  *  1  C(ij+  1 

j=i 


/  £  C 


{£  -  i 
J 


n 

TT  65 

S  ss  ^  X/  c  j  0 


(21) 


If  we  neglect  the  electron  recoil  the  equation  (21) 
is  solved  straightforwardly  and  reads 
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CU}(ke’T) 


=  .  tt  ~=z  exp  |  i  -7  |  |a.(T')|2dT' 

‘s.J-  J=1  /£J  1  2jo  j 


•  exp 


-  ~ | a. ( x) | 2  •  ( a  .  (x ) )  3 
2  J  J 


ct.(t)  =  (-i)Q  .exp  in  .  - 
J  Kj  j  2 


sin(n  . t/ 2 


J 

n  .  /  2 

j 


T-) 


(22) 


Therefore  we  find  that  the  probability  of  finding 
1  .  photons  at  the  time  T  disperses  as  a  Poissonian. 
The  inclusion  of  the  electron  recoil,  even  at  the 
first  order,  is  a  quite  complicated  job  f 1 3]  ,  but 
just  to  give  an  idea  we  have  evaluated  the  average 
number  of  photons  in  the  n-th  mode,  including  the 
first  order  quantum  correction,  thus  finding 


+  °° 


<Z  >  = 
ra 


e  la  (T)l2  +  2  ?  [  ,R  R(t,n.,n  ) 

mm  dn  m  j  jm  e  j  ra 

ra 


-f  |ajW|2  3-^-|aB(t)n|an(t)|  + 


+  2  |£j,m  |a.(x  )|  2  )|am(T)  |  2 


m 


-  2  |  c3.«  I«J<T)|  ^  |a.<x)|>  * 

ra 


sinq  T 
ra 


v  (ni-n?)  I 

“  ejra  “L~  |a.(t)|2(  dk-  . 

n  1  j  '  n  .  j  \  € 

m  j 


(23) 


where  R(T,  7;’7»  )  is  a  somewhat  complicated  function  and  it  will  be 
reported  elsewhere  [l3],  .  Fran  (23)  it  follows  that  in  first  order  in 
£ll  the  laser  modes  evolves  independently  of  each  other  but  becomes 
coupled  for  higher  order  of  the  coupling  parameter.  This  result  should 
be  campared  to  that  of  Ref. [12] where  it  was  shown  that  the  gain  of  the 
FEL  in  the  longitudinal  multimode  hypothesis  is  linked  to  the  Fourier- 
transform  of  the  longitudinal  electron-beam  distribution.  Finally,  we 
want  to  point  out  that  the  above  outlined  multimode  quantum  theory  may 
be  the  basis  for  the  understanding  of  the  growth  of  the  FEL  signal  from 
the  noise  [13]- 
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Table  I  List  of  the  symbols  used  throughout  the  text 


Electron  mass 
Light  velocity 
Planck  constant 


r0  =  e  2 /me  2 

P°,  z 

bL,U’  bL,U 


k  “~k 
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L,  U 


At 

V 


wo 
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a 
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a 

R 

a0  =  !<*o 


-iut 


e 


k 

s 

k 
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Classical  electron  radius 
Electron  longitudinal  mo- 
raentun  and  coordinate 
Annihilation  and  crea¬ 
tion  operators  for  laser 
and  undulator 
Laser  and  undulator  wave  — 
-vectors  in  the  moving 
frame 

Laser  (undulator)  fre¬ 
quency 

Number  of  photons  of 
the  laser  and  the  undulator 
Interaction  time 
Interaction  volume 
2ui  At  P0  /me 
2!i  k2  /m  At 
2TTC2ft0/u)V 
QAt 

n 

I  Qt0  !  ^average  number  of  p ho- 
tons  in  the  input  laser  field 
s-th  laser  mode  wave  vector 
electron  wave  vector 
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SJ 


electron  wave  packet  dis¬ 
tribution  in  the  momentum 

representation 
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